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Motivation 


There are many transport codes for space radiation calculations 
of heavy ion transport: 

Deterministic transport (1 -dimensional): 

HZETRN (from NASA Langley Research Center) p as t 

UPROP (from Naval Research Lab) 


Monte Carlo transport (3-dimensional): 
HETC, HETC-HEDS (from LANL/NASA/ORNL/UTK) 
FLUKA(from high energy physics) 

GEANT4 (from high energy physics) 

MCNP and MCNPX (from LANL) 

PHITS (from Japan/Sweden) 


slower, but better treats 
3-d particle transport 


Are these models very different in typical space radiation calculations? 

How different are 1-d deterministic results from 3-d Monte Carlo (MC) results? 
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Main physics in radiation transport codes 

Let us look at a radiation transport equation in 1 -dimension: 


5J k (E,x) 
t dx 


Jk(E,x) | y Jj(E,x) ^ 
A k (E) j A k j(E) 


5[w k (E)J k (E,x)] 

0E 



due to its fragmentation: 
A k (E)=1/(n*c k (E)) 

total inelastic cross section 
of nuclear fragmentation of k 


Ionization energy loss 
w k (E) = -dE/dx(E) 

Gain of k because a 
heavier particle j can 
produce k: 

A kj (E) = 1/(n*c kj (E)) 

Partial fragmentation 
cross section (j^k) 


Fragmentation cross sections & energy loss 

are the key physics in radiation transport codes 
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Radiation transport codes being considered 


Deterministic transport (1 -dimensional): 

HZETRN (from NASA Langley Research Center) 
UPROP (from Naval Research Lab) 


> 1 995 version used here 


Monte Carlo transport (3-dimensional) : 
HETC, HETC-HEDS (from LANL/NASA/ORNL/UTK) 


FLUKA(from high energy physics) - 
3EANT4 (from high energy physics) 


> 


MCNP and MCNPX (from LANL) 
PHITS (from Japan/Sweden) 


very limited results 


To compare the key physics in radiation transport codes, 
we consider the same 

radiation environment input, geometry & material, 
then compare 

dose-depth curves & particle spectra. 
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Space radiation cases being considered 

3 external environments: 

• Oct. 1989 Solar Particle Event (SPE), 

•Jan. 2005 SPE, 

• 1977 solar minimum Galactic Cosmic Rays (GCR) 

2 materials: Aluminum or CH 2 


2 geometries: 

1) Slab geometry: 
a slab material 

under uni-directional irradiation 


2) Spherical geometry: 

a spherical shell under isotropic irradiation, 
shell thickness 10g/cm 2 
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Oct 1989 SPE, after 40g/cm i A1 

FLUKA at 40.5g/cm 2 


Results for slab geometry: SPE 

Oct 1989 SPE on Aluminum Slab 
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FLUKA, Geant4 & HZETRN are consistent in dose, 
some difference in dose equivalent; 

UPROP dose is lower behind shielding 


HZETRN & Monte Carlo show 
reasonable agreement in proton spectra, 
difference in neutrons, especially at low energies; 
UPROP gives a lower proton spectrum 
(Note: UPROP does not treat neutrons) 



Energy (MeV) 


Oct 1989 SPE, after 80g/cnE A1 



Energy (MeV) 
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Results for SPE: 

new(2005) & old(1995) versions of HZETRN 


2005 version of HZETRN includes 
improved neutron transport , 
especially at low energies (<100 MeV) 


Feb 1956 Webber SPE spectra after 20g/cnT Aluminum 
OLTARIS 2009 » HZETRN 2005 
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Neutron transport in HZETRN 
is still evolving with time 
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Proton transport remains -same 


OLTARIS: https://oltaris. larc. nasa.gov 
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• • HZETRN 1 995 (after Al) 

O- - -O From OLTARIS (after Al) 

HZETRN 2005 (after Al before 30g/cm 2 water) 

From OLTARIS (after Al before 30g/cm 2 water) 
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Energy (MeV) 

HZETRN 2005 reference: 

Heinbockel et al., NASA-TP-2009-21 5560 
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1977 GCR after lOg/crrf Aluminum 


Results for slab geometry: GCR 

1977 GCR on Aluminum Slab 



Geant4 & HZETRN are consistent in dose; 
UPROP dose is lower behind shielding 


HZETRN & Monte Carlo show 
rough agreement in proton spectra, 
difference in neutrons, especially at low energies; 
UPROP proton spectrum is much lower 
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Differential fluence rate in #/(cm / -sr-yr-MeV/u) 


Results for slab geometry: GCR 


1977 GCR after lOg/cm" Aluminum 1977 GCR after 100g/cm~ Aluminum 




HZETRN, UPROP & Geant4 

show reasonable agreements in Oxygen & Iron spectra, 

— > fragmentation cross sections are similar in these models (for O & Fe at least) 
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Results for spherical geometry 



An Aluminum spherical shell in Oct 1989 SPE 


Inner radius 150cm, shell thickness 10g/cm“ 
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An Aluminum spherical shell in 1977 GCR 


Inner radius 150cm, shell thickness lOg/cm 
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Good agreement for the SPE environment 


Agrees with earlier findings based 1 -dimensional transport: 

Lin, Baal I a & Townsend, Radiation Measurements 44 (2009) 

“Variation of space radiation exposure inside spherical and hemispherical geometries” 

• lowest radiation exposure is at the inside wall, 

• highest exposure is at the center of the spherical shell; 

• exposure decreases by a large factor in SPE environments 
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Summary 


We have compared typical space radiation calculations from 
two 1 -dimensional deterministic codes (HZETRN, UROP) 

& two 3-dimensional Monte Carlo codes (FLUKA, Geant4) 

Monte Carlo codes (FLUKA and Geant4) results are mostly consistent, 
HZETRN results are close to Monte Carlo results, except for neutrons 

UPROP results are often quite different from the other 3 codes, 
suggesting the need of improvements (e.g. by treating neutrons) 

Radiation exposure at different locations inside a spherical shell: 

4 models give consistent results, 

earlier result of Lin, Baalla & Townsend is confirmed 

To identify the exact physics causing the differences in the model results 
is very useful but will require more efforts 
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